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Abstract 

Bose-Einstein correlations between like-charge pions are studied in hadronic final 
states produced by e + e~ annihilations at center-of-mass energies of 172 and 183 
GeV. Three event samples are studied, each dominated by one of the processes 
W + W~ — > qqfi^v, W + W _ — ► qqqq, or (Z /j)* — ► qq. After demonstrating the 
existence of Bose-Einstein correlations in W decays, an attempt is made to determine 
Bose-Einstein correlations for pions originating from the same W boson and from 
different W bosons, as well as for pions from (TP/j)* — ► qq events. The following 
results are obtained for the individual chaoticity parameters A, assuming a common 
source radius R: 

A samo = o.63 ±0.19 ±0.14, 
A difr = 0.22 ±0.53 ±0.14, 
A z * = 0.47 ±0.11 ±0.08, 
R = 0.92 ± 0.09 ± 0.09 fm. 

In each case, the first error is statistical and the second is systematic. At the current 
level of statistical precision it is not established whether Bose-Einstein correlations 
between pions from different W bosons exist or not. 
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1 Introduction 



In reactions leading to hadronic final states Bose-Einstein correlations (BEC) between 
identical bosons are well known. These correlations lead to an enhancement of the number 
of identical bosons over that of non-identical bosons when the two particles are close 
to each other in phase space. Experimentally this effect was first observed for pions 
by Goldhaber et al. For recent reviews see, for example, reference 0. In e + e~ 
annihilations at center-of-mass energies of 91 GeV, BEC have been observed for charged 
pion pairs for K°K° pairs JTHKj and also for K ± K ± |T|. 



In the present paper we report on an investigation of BEC for charged pions between 
e + e~ reactions at center-of-mass energies of 172 and 183 GeV, above the threshold for W- 
pair production. The analysis is motivated by the question of whether BEC for pions from 
different W bosons exist or not. Theoretically this question is still not settled fl~2|,[T3" 



However, if such correlations do exist, this could bias significantly the measurement of the 
W boson mass in fully hadronic W-pair events [P^, P^TB[|. The DELPHI collaboration 
has published a measurement, at y/s = 172 GeV, of BEC between pions originating from 
two different W bosons ||17|| , in which basically BEC in W + W _ — > qqiVg events were 



subtracted from those of W + W _ — > qqqq events. The aim of the present analysis is to 
analyse BEC for fully hadronic W-pair events (W + W~ — > qqqq), semileptonic W-pair 
events (W + W~ — > qqivg), as well as non-radiative (Z°/j)* —>■ qq events. After having 
established BEC in hadronic W decays, BEC are investigated separately for three classes 
of pions: those originating from the same W boson, those from different W bosons and 
those from non-radiative (Z° /'-/)* — > qq events. Note that in this analysis, tracks are not 
assigned to jets or W-bosons and no kinematic fits are needed. 



BEC between identical bosons can be formally expressed in terms of the normalised 



function 

C(Q) = p2 ( pi i p2 "> = ^ d2(T I ( da da 
Pi(Pi)pi(P2) dp 1 dp 2 / \dp 1 dp 2 

where a is the total boson production cross section, piQoj) and da/dpi are the single- 
boson density in momentum space and the inclusive cross section, respectively. Similarly 
P2(Pi,P2) an d d 2 a/dpidp2 are respectively the density of the two-boson system and its 
inclusive cross section. The product of the independent one-particle densities Pi(pi)pi(p2) 
is referred to as the reference density distribution, to which the measured two-particle 
distribution is compared. The inclusive two-boson density ^2(^1,^2) can be written as: 

p2(puP2) = pi(pi)pi(p2) + K 2 (pi,p 2 ) , (2) 

where K 2 (pi,p 2 ) represents the two-body correlations. In the simple case of two identical 
bosons the normalised density function C(Q), defined in Eq. [IJ, describes the two-body 
correlations. Thus one has 

C(Q) = I+K2 (PUP2) , (3) 

where K2 {Pi,P2) = K 2 (pi, p 2 ) / [pi{pi) Pi{p2)} is the normalised two-body correlation term. 
Since BEC are present when the bosons are close to one another in phase space, a natural 
choice is to study them as a function of the Lorentz invariant variable Q defined by 

Q 2 = -{pi-p2? = Mi-^ 2 , 

which approaches zero as the identical bosons move closer in phase space. Here p% is 
the four-momentum vector of the zth particle, \x is the boson mass (here m,) and Mf 
is the invariant mass squared of the two-boson system. Ideally the reference sample 
should contain all correlations present in the sample used to measure p{pi,p 2 ), other than 
the BEC , such as those due to energy, momentum and charge conservation, resonance 
decays and global event properties. In this analysis, the reference is chosen to be a 
sample of unlike-charge pairs of pions from the same event. Since the presence of the 
resonances u, Kg, rj, rj', p°, fo and f 2 in the unlike-charge reference sample leads to 
kinematic correlations which are not present in the like-charge sample, the unlike-charge 
sample has to be corrected for this effect using simulated events. 

Assuming a spherically symmetric pion source with a Gaussian radial distribution, the 
correlation function C(Q) can be parametrised [|TJ by 

C(Q) = N(1 + U(Q)\e- Q2R2 ) (1 + 5Q + eQ 2 ), (4) 

where R is the radius of the source and A represents the strength of the correlation, 
with < A < 1. A value of A = 1 corresponds to a fully chaotic source, while A = 
corresponds to a completely coherent source without any BEC. The function fn{Q) is 
the probability that a selected track pair is really a pair of pions, as a function of Q. 
The additional empirical term (1 + 5Q + eQ 2 ) takes into account the behaviour of the 
correlation function at high Q values due to long-range particle correlations (e.g. charge 
and energy conservation, phase-space constraints), and iV is a normalisation factor. 

The structure of the paper is as follows. Section 2 contains a brief overview of the 
OPAL detector, the event and track selections as well as Monte Carlo models. In section 3 
the analysis of the data is described. BEC are investigated for (ZP/j)* — > qq, W + W~ — ► 




qqiVi and W + W~ — > qqqq events. After establishing BEC in hadronic W-events the 
chaoticity parameter for BEC between the decay products from the same W, A same , and 
from different W bosons, A dlff , are determined. Finally, section 4 summarises the results 
obtained. 



2 Experimental Details 
2.1 The OPAL detector 

A detailed description of the OPAL detector has been presented elsewhere [ITS) and there- 



fore only the features relevant to this analysis are summarised here. Charged particle 
trajectories are reconstructed using the cylindrical central tracking detectors which con- 
sist of a silicon microvertex detector, a high-precision gas vertex detector, a large- volume 
gas jet chamber and thin z-chambers Q. The entire central detector is contained within a 
solenoid that provides an axial magnetic field of 0.435 T. The silicon microvertex detector 
consists of two layers of silicon strip detectors, allowing to measure at least one hit per 
charged track in the angular region | cos#| < 0.93. It is surrounded by the vertex drift 
chamber, followed by the jet chamber, about 400 cm in length and 185 cm in radius, that 
provides up to 159 space points per track and also measures the ionisation energy loss of 
charged particles, dE/dx. With at least 130 charge samples along a track, a resolution 
of 3.8% is achieved for the dE/dx for minimum ionising pions in jets [|l^, 2(J. The z- 



chambers, which considerably improve the measurement of charged tracks in 6, follow the 
jet chamber at large radius. The combination of these chambers leads to a momentum 
resolution of cr p /p 2 = 1.25 x 10~ 3 (GeV/c) _1 . Track finding is nearly 100% efficient within 
the angular region | cos^l < 0.92 . The mass resolution for Kg — > 7r + 7r~, related to the 
resolution in the correlation variable Q, is found to be a = 7.0 ± 0.1 MeV/c 2 M. 



2.2 Data selection 



This study is carried out using data at e + e~ center-of-mass energies of 172 GeV and 
183 GeV with integrated luminosities of approximately 10 pb _1 and 57 pb _1 , respec- 
tively. Three mutually exclusive event samples are selected: a) the fully hadronic event 
sample, W + W~ — > qqqq, where both W bosons decay hadronically; b)the semileptonic 
event sample, W + W _ — > qqiVi, where one W decays hadronically and the other decays 
semileptonically ; and c) hadronic non-W events (ZP/j)* — > qq, referred to here as the the 
non-radiative (Z°/j)* event sample in this analysis. Throughout this paper, a reference 
to W + or its decay products implicitly includes the charge conjugate states. 



1 The OPAL right-handed coordinate system is defined such that the origin is at the geometric centre 
of the jet chamber, z is parallel to, and has positive sense along, the c~ beam direction, r is the coordinate 
normal to 2, 8 is the polar angle with respect to +2 and <fi is the azimuthal angle around 2. 



2.2.1 Selection of the fully hadronic event sample W+W 



qqqq 



The selection of fully hadronic W + W~ — > qqqq events is performed in two stages using 
a preselection based on cuts followed by a likelihood-based selection procedure. Fully 
hadronic decays, W + W~ — > qqqq are characterised by four or more energetic hadronic jets 
and little missing energy. A preselection using kinematic variables removes background 
predominantly from radiative (ZP/j)* — > qq events. Events satisfying the preselection 
criteria are subjected to a likelihood selection, which discriminates between signal and 
the remaining four-jet-like QCD background. 

At 172 GeV, several variables based on the characteristic four-jet-like nature, momen- 
tum balance and jet angular structure, are used to distinguish W + W~ — > qqqq events 
from the remaining background and to construct the likelihood. The details of the selec- 
tion at 172 GeV are described in appendix B of [^1]]. The signal and background situation 
at 183 GeV is similar to the one at 172 GeV. For this reason, no new selection strategy 
was developed and the event selection at 183 GeV is just a reoptimised version of the 



selection at 172 GeV. The details of the selection at 183 GeV are described in |[22|| . At 
183 GeV, no cut was applied against Z°Z°events. 

Overall, there is a background of 11.6% from (Z°/7)* — > qq events and a contribution 
of 2.1% from e + e~ — > Z°Z° — > qqqq events. No selection for W + W~ — > qqqq events is 
applied to events selected as W + W~ — > qq£&£ events. 



2.2.2 Selection of the semileptonic event sample W+W — ► qq£is£ 

W + W~ — ► qqez7 e and W + W~ — > qqu^V events are characterised by two well-separated 
hadronic jets, a high-momentum lepton and missing momentum due to the unobserved 
neutrino. In W + W _ — > qqrV T the r lepton gives rise to a low-multiplicity jet consisting 
of one or three tracks. The tracks from r decay are not used in in the BEC studies. Cuts 
are applied to reduce the background from radiative (Z°/j)* — > qq events. A likelihood 
is formed using kinematic variables and characteristics of the lepton candidate to further 
suppress the background from (Z /^)* — > qq events. The details of the selection at 172 
GeV are given in appendix A of [^J . 

The W + W~ — > qq£p£ event selection for the 183 GeV data is a modified version of 
the 172 GeV selection. At 183 GeV, a looser set of preselection cuts is used since the 
lepton energy spectrum is broader due to the increased boost and the set of variables 
used in the likelihood selections is modified. In the W + W _ — > qqrV T sample there 
is a significant background from hadronic decays of single W events (e + e~ — > WeU e ) 
and an additional likelihood selection is used to reduce this background. This is only 
applied to W + W~ — > qqrz7 T events where the tau is identified as decaying in the single 
prong hadronic channel. Finally, in order to reduce Z°Z° contribution, events passing 
the W + W~ — ► qqez/ e likelihood selection are rejected if there is evidence of a second 
energetic electron. A similar procedure is applied to the W + W _ — > qq/il 7 ^ selection. The 
details of the selection at 183 GeV are given in [p2]| . There is a background of 3.5% from 
(Z°/j)* — > qq events, 1.0% from W + W~ — > qqqq events, 1.3% from single W events and 
0.8% from Z°Z° -> qqii events. 



2.2.3 Selection of the non-radiative event sample (ZP/~y)* 



qq 



Here, an extension of the selection criteria defined in is used, which starts by selecting 
hadronic events defined as in 



25 



To reject background from e + e _ — > t + t~ and 77 — > qq 
and to ensure that the events are well contained in the OPAL detector one requires that 
the event has at least seven charged tracks with transverse momentum p t > 150 MeV/c 
and that the polar angle of the thrust axis lies within the range | cos8t\ < 0.9. To reject 
events with large initial-state radiation, one requires \fs — y/s' < 10 GeV, where \fs i is the 



effective invariant mass of the hadronic system [23]. For the suppression of the W + W 



background one requires that the events are selected neither for the semileptonic nor for 
the fully hadronic W + W~ samples described above. The cut in the relative likelihood for 
vetoing W + W~ — > qqqq events is looser than in the W + W~ — > qqqq event selection. After 
selection, there is a residual background of 3.8% from W-pair events and a contribution 
of 0.3% from Z°Z° events. 



2.2.4 Pion selection and Event samples 



Note that the three event selections result in completely independent event samples with- 
out any overlap. After the event selection the following cuts are applied to all tracks, for all 
three event samples. A track is required to have a transverse momentum p t > 0.15 GeV/c, 
momentum p < 10 GeV/c and a corresponding error of a p < 0.1 GeV/c. Only tracks with 
polar angles 9 satisfying | cos^l < 0.94 are considered. The probability for a track to be a 
pion is enhanced by requiring that the pion probability P^ from the dE/dx measurement 
is Y n > 0.02. Pion-pairs from a Kg decay are rejected using the Kg finder described in ||. 
This algorithm rejects 31% of the unlike-charge pion pairs coming from a Kg decay. Since 
less than 11% of the rejected pairs do not originate from a Kg, this cut does not intro- 
duce a significant bias in the Q-distribution. Finally, events with fewer than five charged 
selected tracks are rejected. The number of events retained, as well as the number of 
background events evaluated from Monte Carlo simulation is given in table |T|, for all three 
event samples. 



event sample 


number of selected events 


expected background events 


172 GeV 


183 GeV 


172 GeV 


183 GeV 


W + W _ — > qqqq 


55 


327 


9.5 ±0.5 


43.6 ±2.4 


W+W" -> qq£u e 


45 


326 


2.1 ±0.5 


23.1 ±2.4 


(Z°/ 7 )* 


214 


1009 


8.1 ± 1.7 


43.2 ±4.9 



Table 1: Number of retained events and number of background events predicted for the 
three event samples, separately for 172 GeV and 183 GeV. 



2.3 Monte Carlo models 



A number of Monte Carlo models are used to model (Z°/j)* — > qq, W + W~ — > qqiVg or 
W + W _ — > qqqq events. For the W + W _ — > qqqq event sample the simulated events are 
also used to determine the fraction of track pairs coming from the same or different W 



bosons. The Monte Carlo samples are generated at e + e~ center-of-mass energies of 172 
and 183 GeV in proportion to the corresponding integrated luminosities. The production 

Non-radiative decays (Z /^)* — > qq as well as 



Koralw uses the same string 



of W-pairs is simulated using Koralw 
the Z°Z° and We77 e events are simulated with Pythia 
model as Pythia for hadronisation. For systematic error studies the event generator 
Herwig |23| , which employs a cluster hadronisation model, is also used. All these Monte 
Carlo samples discussed above are generated without BEC. In addition W-pair events are 
also simulated with BEC included [I5|, using Pythia Q The algorithm introduces BEC 
via a shift of final-state momenta among identical bosons. For these events two samples 
are generated: In the first, BEC are simulated for all pions in the event, both from the 
same and from different W bosons. In the second sample, BEC are simulated only for 
pions originating from the same W boson. 



3 Analysis 

Using the tracks that pass the selection of section |2.2.4|, the Q-distributions are determined 



for like-charge pairs as well as for unlike-charge pairs. The correlation function C(Q) is 
then obtained as the ratio of these Q-distributions. Coulomb interactions between charged 
particles affect like- and unlike-charge pairs in opposite ways and modify the correlation 
function. We therefore apply the following correction to the correlation function, 

C oarr (Q) = x(Q)C wcmr (Q), (5) 

where 

x(Q) = t—=^> (6) 

and where rj = a m^/Q with a the fine-structure constant and the mass of the charged 
pion P5 . The Coulomb correction factor x{Q) is about 17% in the first Q bin, 5% in the 



second bin and 1% in the tenth bin, with a bin size of 0.08 GeV/c 2 (see Fig. [I] for the 
definition of the bins). The Monte Carlo simulations do not contain Coulomb effects, so 
the Monte Carlo distributions are not corrected by Eq. |5|. 

Structure in the unlike-charge samples due to resonance production is corrected using 
Monte Carlo. For this, the Q-distribution is obtained for unlike-charge pair combinations 
taken exclusively from the decay products of Kg mesons and the resonances u, rj, rf, 
p°, fo and f2 as produced in the Monte Carlo. The production of resonances has only 
been measured at e + e~ center-of-mass energies around the Z° peak and not at energies 
above the Z° peak. Jetset |3(| describes the production of resonances around the Z° 



peak quite well, although not perfectly in all cases [|3T]. To estimate the contribution 



for each resonance to the Q-distribution at LEP 2 energies, the Q-distribution for each 



resonance is multiplied by the ratio of the measured production rate at LEP [31] and 
the corresponding rate in Jetset. The main contributions come from Kg, uj, p° and 
rj mesons. The Q-distribution for the resonances, thus obtained, is then scaled to the 

2 The model parameters controlling BEC in Pythia are taken to be MSTJ(51)=2, MSTJ(54)=-1, 
MSTJ(57)=1, PARJ(92)=1.0, PARJ(93)=0.4, MSTJ(52)=9, PARJ(94)=0.275, and PARJ(95)=0.0, as 
suggested by the authors of [Hq] . 



number of selected events and subtracted from the experimental unlike- charge reference 
Q-distribution. These corrections are made for each event selection separately. They are 
typically 5 — 10% for small Q, falling rapidly for Q > 0.8 GeV/c 2 . The three unlike-charge 
distributions, before the correction, and the expected signal from resonance decays are 
shown in Fig. [I]. 

The resulting experimental correlations C(Q) are shown, for the three event samples 
separately, in Fig. |[ The data in all three distributions exhibit a clear enhancement at 
low Q, consistent with the presence of BEC. 



3.1 Fit to establish BEC in W-pair events 

The measured distributions cannot be directly compared with the parametrisation of 
Eq. |] since, in general, each distribution has contributions from several physical pro- 
cesses that may have different BEC. To illustrate the situation, consider the hadronic 
W-pair events. They have as their main contribution (see Table [3]) the correlations from 
pions coming from hadronic W decays. They contain, however, also contributions from 
background events, i.e. (Z°/j)* — > qq events. Thus 

u W) jyw_w + N z*_ ' V) 

where and N± ± are the numbers of like-charge track pairs for the class of pions 

from W + W~ — > qqqq events and for the class of pions from the background sample of 
(Z°/j)* —>■ qq events. The variables jVjX^ and N+_ are defined analogously for unlike- 
charge pairs. Eq. |^ can be rewritten as 

C had (g) = pWW (Q) C ^ {Q) + pWW {Q)) C g d( g )i (8) 



where C q ™((5) and C^* d (Q) are the BEC for the class of pions from W + W" — ► qqqq 
events and for the class of pions from the sample of (Z° /'-/)* — ► qq events, the main 
background in the hadronic selection. P^^(Q) = /(N^^ + iV?_) is the fraction of 
unlike-charge pion pairs at a given Q which originate from a W-pair event in the hadronic 
event sample. Here and in the following, the small number of Z°Z°events are not counted 
as background but as signal, since their properties with regard to BEC should be quite 
similar. 

The experimentally determined correlations for the other two event samples can be 
written as: 

c scmi (Q) = pZM c*{Q) + (i - pZAQ)) c z \Q), (9) 

for the W + W~ — > qqiVi event sample and 

C" d -(Q) = P n z o ;_ rad (Q) C Z *(Q) + (1 - P n z * n _ rad (Q)) C*«(g) (10) 

for the non-radiative (Z°/j)* event sample. The notation in these equations is analogous 
to that of Eq. [7|. C semi (Q) and C z * (Q) are the BEC for the two pion classes from 
W + W~ — > qqiVi and non-radiative (Z° /'-/)* —>■ qq events, respectively. The definition of 
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Figure 1: The points show the unlike-charge pion pair distribution, before the correction 
for resonance production, for the three different event selections; a) for the fully hadronic, 
b) for the semileptonic, and c) for the non-radiative event selection. The filled histogram 
is the expected contribution from resonances, where both tracks of a pion pair come from 
the same resonance. 
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Figure 2: The correlation function for like-charge pairs relative to unlike-charge pairs for 
three event selections; a) C ha,d (Q) for the fully hadronic, b) C semi (Q) for the semileptonic, 
and c) C non ~ rad ' (Q) for the non-radiative event selection. The Coulomb-corrected data 
are shown as solid points together with statistical errors. The curves are the result of the 



simultaneous fit discussed in Sect. 3.2. 



the relative fractions P^^iQ), P™ m i(Q) an d Pnon-nd(Q) * s gi ven i n table 0. They are 
taken from a Monte Carlo simulation which does not contain BEC as discussed in section 



273]. These probabilities are global properties of the events and depend little on whether 
BEC are assumed or not. The small number of single- W events in the semileptonic event 
sample are treated as signal events. 

The hadronic W-pair sample contains a sizeable number of (Z /^)* background events. 
Due to the selection cuts suppressing (Z°/j)* events in the hadronic W-pair sample, the 
remaining (Z°/7)* events have different event shapes and multiplicities from those in 
the main non-radiative (Z°/j)* event sample. Since BEC depend on event shape and 
multiplicity {§], the correlation function for (Z°/j)* events selected as hadronic W-pairs, 
Chad (Q) > i s expected to be different from that for the main non-radiative selection, C z * (Q). 
To take these differences into account, the parameters A and R in the correlation functions 
C^ad(Q) and C z * (Q) are not taken to be equal but those in C z ^ d (Q) are adjusted according 
to the different event topology. In order to estimate this correction, the W + W~ — > 
qqqq selection described in |2]2|, which contains no direct center-of-mass energy dependent 
variables, is applied to data taken at LEP. A simultaneous BEC fit is applied to both 
events selected as W + W~ — > qqqq events and events which are not selected as W + W~ — ► 
qqqq events. The differences obtained in A and R are used here f\ to take differences 
in the correlation function C^ d (Q) and C z * (Q) into account. Due to high purity of 
the semileptonic and non-radiative (Z°/j)* selections, no adjustment is applied to the 
correlation functions of events selected as background in C semi (Q) and C non ~ rad '(Q). For 
W + W _ — > qqqq events selected as fully hadronic events and W + W~ — > qqqq events 
selected as non-radiative (Z°/7)* events the same BEC are assumed. The effect of this 
assumption will be described with the systematic errors. 
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Figure 3: The purities a) P^^(Q), b) P^ mi (Q) and c) P^^ iad (Q) as obtained from 
Monte Carlo simulations. 



3 For the function C^ ad (Q) the absolute A value is reduced by 0.094 and the absolute R value is 
increased by 0.097 fm relative to the corresponding parameters of C z (Q), with A z kept as a free 
parameter in the main BEC fit. 



Probability definition 


Prob. that H — track pair 


pWWfpN _ N^(Q) 
Had (Q)+N*_(Q) 


originates from W + W~ — > qqqq process, 
in the hadronic event selection. 




originates from (Z /^)* — > qq process, 
in the hadronic event selection. 


P w (n] _ N ?-(Q) 

beiru v -v / N+_(Q)+lv + _(Q) 


originates from W + W~ — > qqiVg process, 
in the semileptonic event selection. 


P z* rn N _ 

non-radV"W jyww (Q)+jV~_ (Q) 


originates from (Z°/j)* — > qq process, 
in the non-radiative event selection. 


pWW _ i _ pZ* /n\ 


originates from W + W~ — > qqqq process, 
in the non-radiative event selection. 


jysamc W/zn'i 

^had — ArpmcW ( Q) +A rdiffw ( Q) +7V z*_ ( Q) 


originates from the same W, 
in the hadronic event selection. 


jusamc W//n\ 

psamc ((O") ^ — 

-'non-radWJ — JV«? ieW (Q)+jyf"? w (Q)+iV|!.(Q) 


originates from the same W, 
in the non-radiative event selection. 



Table 2: Definition and meaning of the various probabilities concerning unlike-charge 
track pairs, used in Eqs. |8| - [10] and [12| - |14] and illustrated in Figs. |3| and [|. 



The unknown correlation functions C q ™(<5), C°^(Q) and C z (Q) are parametrised 
using Eq. |j. The parameters are determined in a simultaneous fit to the three experimental 
distributions shown in Fig. p.2| . A common source radius R is used for all event classes, 
while the parameter A is allowed to be different. There is justified by the typical separation 
between the W + and W~ decay vertices is smaller than 0.1 fm at LEP 2 energies, much 
smaller than the typical hadronic source radius of R ~ 1 fm \I2\, justifying equal source 
radii for the W + W~ — > qqqq and the W + W~ — ► qq^ event classes. The source radius 
for e + e~ annihilations into hadrons has been measured up to 90 GeV and no evidence has 
been found for an energy dependence ||. For this reason R is assumed to be the same 
at higher energies and the same source radius is also used for the (Z°/j)* —>■ qq event 
class. Separate fits to the distributions show also consistent radii for the different event 
selections. The pion probability fn{Q) is taken from Monte Carlo. At small values of Q it 
is about constant ~ 0.84 and varies only weakly with Q for all channels. The long-range 
parameters are expected to be different for the W + W~ — > qqqq, W + W~ — > qqiVi and 
(Z°/7)* — > qq event class, due to kinematic and topological differences. The results for the 
thirteen free parameters in the fit are given in Table |3|. The fit is made in the full range of 
0.0 < Q < 2.0 GeV/c 2 . In the distributions of Fig. [2] the same particles contribute many 
times, in different bins of Q, which introduces bin-to-bin correlations. These are taken 
into account in the fit. All three experimental distributions are well described by the fit, 
the xVd.o.f. is 76.1/62. 



Parameter 


W + W — > qqqq 


W+W" -> qq£u E 


(Z°/7)* 


R (fm) 




0.91 ±0.11 ±0.10 




A 


0.43 ±0.15 ±0.09 


0.75 ±0.26 ±0.18 


0.49 ±0.11 ±0.08 


N 


0.86 ±0.04 ±0.04 


0.79 ±0.08 ±0.08 


0.86 ±0.05 ±0.04 


5 


0.12 ±0.10 ±0.10 


0.29 ±0.23 ±0.24 


0.13 ±0.11 ±0.08 


e 


-0.04 ±0.05 ±0.06 


-0.09 ±0.10 ±0.11 


-0.02 ±0.05 ±0.04 



Table 3: Result of the simultaneous fit. The first error corresponds to the statistical 
uncertainty the second to systematics. 



3.2 Fit to establish BEC in same and different W bosons. 

In this section BEC are investigated separately for pions originating from the same W 
boson and for pions from different W bosons. The correlations for the fully hadronic event 
sample (Eq. |7|) are written as 

ArsameW _j_ /vydifTW _j_ AfZ* 

jV^ameW _|_ jydiflfW _|_ jyZ* ' \ lL ) 

where A^ a ™ eW , iV^ w and N^ ± are the number of like-charge track pairs for the class of 
pions from the same W boson, different W bosons and from (7P/j)* — > qq events. The 
variables A^^. a ™ eW , jy^z and iV?_ are defined, in a similar way, for unlike-charge pairs. 
Eq. [11] can be rewritten as 

c had (Q) = P£T(Q) c samc (Q) + pZAQ) cEd(Q) 

±(i - PITiQ) - PgiiQ)) c diS (Q), (12) 



where C same (Q),C dlS (Q) and C z (Q) are the BEC for the class of pions from the same 
W boson, different W bosons and from (ZP/j)* — > qq events. The variables Pll™ e (Q) and 
P^adiQ) are defined in Table || Likewise, the experimentally determined correlations for 
the other two event samples can be written as: 

C^(Q) = P s W mi (Q) C"(Q) + (1 - P^M) C z \Q) (13) 

for the W + W~ — > qq£z/£ event sample and 

C— d '(Q) = P n s a n m ° ra d(Q) C samc (Q) + C Z *(Q) 

+(i - p^UQ) - pL-UQ)) c diS (Q) (14) 

for the non-radiative Z* event sample. The definition of the variables P^ S (Q) , Pha,d(Q)> 
PZiiQ), P^n-UQ) , and P n z o ;„ rad (Q) is also given in table g 

By simultaneously fitting Eq's. 0, and [TJ] to the experimental distributions in 
Fig. |, the BEC for the three pion classes C same (Q), C diS (Q) and C z * (Q) are determined. 
Again, the probabilities P n s T-rad(Q), i^-™d(Q), ^(Q) and ?1(Q) are 

taken from Monte Carlo simulations not containing BEC, as discussed in section [2.3|. 



The functions Plon-rad(Q) and P^ e (Q) are shown in Fig. Their properties contain 
only information from unlike-charge pion pairs and are therefore independent of BEC. The 
effect of possible variations of the function P^^ e {Q), if BEC are assumed in the Monte 
Carlo, is discussed in section |3]3|. The unknown correlation functions C samc (Q), C dlS (Q) 
and C z * (Q) are parametrised using Eq. f| with different A for the three event classes. As 
before a common source radius R is used for all event classes. For the correlation function 
ChldiQ) the parameters A and R are adjusted like in section 3A. Based on Monte Carlo 



studies the long range parameters 5 dlS and e for the correlation function C dlS (Q) are 
taken to be zero. This is equivalent to the assumption that colour reconnection effects 
do not influence the Q distributions. The free fit parameters are then determined in a 
simultaneous fit to the three experimental distributions shown in Fig. [| The results for 
the eleven free parameters in the fit are given in Table [|. The fit is made in the full 
range of 0.0 < Q < 2.0 GeV/c 2 . The fit result is given in Fig. ^. All three experimental 
distributions are well described by the fit (x 2 /d.o.f. is 76.4/64). The correlation between 
the parameters A dlff and A same , with a coefficient of —0.52, is shown in Fig. |^. 



Parameter 


same W 


diff W 


(Z°/ 7 )* 


R (fm) 




0.92 ±0.09 ±0.09 




A 


0.63 ±0.19 ±0.14 


0.22 ±0.53 ±0.14 


0.47 ±0.11 ±0.08 


N 


0.83 ±0.05 ±0.07 


1.00 ±0.01 ±0.00 


0.87 ±0.04 ±0.04 


5 


0.21 ±0.15 ±0.19 


zero assumed 


0.11 ±0.11 ±0.07 


e 


-0.07 ±0.07 ±0.08 


zero assumed 


-0.01 ±0.05 ±0.02 



Table 4: Result of the simultaneous fit distinguishing pions from the same and different 
W boson. The first error corresponds to the statistical uncertainty, the second one to 
systematics. 
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Figure 4: The probability that both tracks of an unlike-charge track pair in a fully 
hadronic event originate from the same W boson, -Pno™-rad(Q) (upper plot) and Pl^°(Q) 
(lower plot), as obtained from Monte Carlo simulations. The histogram is the result for 
the case that no BEC are assumed. The dashed and dotted histograms are the results 
for the case that BEC are simulated for all pions or only for pions originating from the 
same W boson, respectively. 




Figure 5: Correlation between A dlfr and A samc . The contour shows the 67% confidence 
level. The best value obtained in the fit is given by the cross. The lines for A dlff = and 
A diff = A samc are also indicated. 



3.3 Systematic Errors 



The following variations in the analysis are considered to obtain the systematic error, 
which affect the fit result in both fit methods. The systematic errors are listed in Table [| 
and [| together with their quadratic sums to give the final systematic error. 



1. Variation of the resonance production. In the main analysis, the distortion of the 
unlike-charge pairs due to resonances was taken into account by subtracting the 
resonance Q distribution from unlike-charge pair Q distribution. This method is 
based exclusively on Z° data, since no measurements of resonance production at 
LEP 2 are available. Thus for systematics the correction factors are varied within 
two standard deviations of the experimental resonance production cross sections. 
The maximum differences in the fit for each resonance were added in quadrature. 
Several variations were made for this systematic check, therefore no x 2 is given. All 
fits are of good quality. 

2. Double-hit resolution. Unlike-charge pairs are bent in a magnetic field in opposite 
directions, whereas like-charge pairs are bent in the same direction. Therefore like- 
charge pairs at very low Q are less well reconstructed. Monte Carlo studies indicate 
the presence of such effects for pairs with a Q less than 0.05 GeV/c 2 . For systematics, 
the fit is repeated in the range between 0.05 < Q < 2.0 GeV/c 2 . 

3. Use of the Herwig Monte Carlo. To determine the purities and the correction for 
resonances in the unlike-charge sample the Herwig Monte Carlo was used. 

4. The probability functions are obtained from Monte Carlo simulations where BEC 
are simulated |16| for all pions, both from the same and from different W bosons. 



5. The probability functions are obtained from Monte Carlo simulations where BEC 
are simulated |16| only for pions from the same W. 



6. Long-range correlations. The fit is repeated with e = 0. 

7. Different topology of the (Z°/j)* — > qq background in fully hadronic selected events. 
The difference of the (Z°/j)* —>■ qq events in the hadronic and non-radiative Z* 



samples is taken into account in the main analysis (see section |3.2| ). The events 
selected at LEP energies as W + W~ — > qqqq events and as non-radiative events, 
which are used for this correction are statistically limited. Therefore the parameters 
governing the correction factor for the correlation function C^ d (Q) are varied within 
their statistical error (A ± 0.04 and R± 0.057 fm) and the largest deviation is taken 
as the systematic error. 



In addition, the effect on uncertainties from arising the knowledge of cross-section is 
examined. The cross-sections for W-pair production processes as well as the cross-section 
for non-radiative (Z /^)* processes are varied within their experimental uncertainties. 
The impact on the final result is negligible. Furthermore, differences between W + W~ — > 
qqqq events selected as hadronic events and selected as non-radiative events are also 
considered. These variations introduce only small changes in the results. 
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4 


+0.01 


-0.02 


-0.02 


-0.01 


75.7/62 


5 


< 0.01 


-0.02 


-0.03 


< 0.01 


75.9/62 


6 


+0.07 


-0.03 


-0.13 


-0.02 


78.2/65 


7 


< 0.01 


< 0.01 


-0.01 


< 0.01 


76.3/62 


total 


0.10 


0.09 


0.18 


0.08 





Table 5: The effect of the systematic variations studied (discussed in Sect. p.3| ) on the 
variables R, A w+W ^ q ™, _\w+w ^qqeu e anc j ^z*_ rp^g j ag ^ co ^ umn ghowg the quality of 

the corresponding fit. 





# (fm) 


^same 


A diti 


\ z 


rVd.o.f. 


Reference 


0.92 + 0.09 


0.63 + 0.19 


0.22 + 0.53 


0.47 + 0.11 


76.4/64 


Variation 


5R (fm) 


^same 


6X m 






1 


±0.07 


±0.09 


±0.07 


±0.07 




2 


< 0.01 


-0.05 


< 0.01 


+0.03 


74.4/64 


3 


+0.01 


+0.04 


+0.03 


< 0.01 


94.5/64 


4 


+0.01 


-0.02 


-0.10 


< 0.01 


76.0/64 


5 


< 0.00 


-0.04 


-0.03 


< 0.01 


76.2/64 


6 


+0.05 


-0.08 


< 0.01 


-0.01 


77.6/66 


7 


< 0.01 


< 0.01 


-0.05 


< 0.01 


76.5/64 


Total 


0.09 


0.14 


0.14 


0.08 





Table 6: The effect of the systematic variations studied (discussed in Sect. |3.3| ) on the 
variables R, A same , A dlfT and A z *. The last column shows the quality of the corresponding 
fit. 



3.4 Q-based separation of BEC contributions 



The experimental BEC for pure classes of a) tracks from different W bosons C dlS (Q) , b) 
tracks from the same W boson C same (Q), as well as c) tracks from (Z /^)* — *> qq events 
C z * (Q) can be obtained directly from Eq. ^ - TA by solving the equations for these three 



unknown functions for each bin of Q, using the fractions from Table |2|. The resulting 
distributions are shown in Fig. |5|. A comparison of data and Monte Carlo without BEC 
shows that there is a clear signal at small Q for pions originating from (Z /^)* — *> qq events 
(Fig. ^ c). The data for pions from the same W boson show a larger enhancement than 
the corresponding simulation (Fig. b). At the current level of precision, it cannot be 
established whether BEC between pions from different W bosons exists or not (Fig. |6| a), 



in agreement with the result of the simultaneous fit of Sect. 3.2. Note that the errors 



of the three unfolded distributions are highly correlated with each other by construction. 
For Q values larger than about 0.4 GeV/c 2 , the distribution for C dlS (Q) is consistent with 
being constant in both MC and data. 
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Figure 6: Correlation functions for the unfolded classes. The data points show the 
experimental distributions for a pure sample of a) pions originating from different W 
bosons C dlS (Q), b) pions originating from the same W boson C same (Q) and c) pions from 
(Z° /'-/)* — > qq events. The errors are the statistical uncertainties and are correlated be- 
tween the three classes. The open histogram in a) is the result, for pions from different 
W bosons, of a simulation including BEC between pions from different W bosons, the 
cross-hatched histogram the corresponding result for a simulation with BEC for pions 
from the same W boson only. The open histogram in b) shows the result, for pions from 
the same W boson, of a simulation including BEC for pions from the same W boson and 
the hatched histogram the corresponding result for no BEC at all. The hatched histogram 
in c) corresponds to a simulation with no BEC at all. 



3.5 Consistency check 



In the analysis described above, the resonances were subtracted using Monte Carlo in- 
formation and long-range correlations were taken into account by the empirical factor 
(1 + SQ + tQ 2 ) in the correlation function of Eq. As an alternative, we study here the 
double ratio 

pjDATA ]yMC 
C (Q) = J^DATA ljqMC~ ( 15 ) 

where the Monte Carlo events are generated without BEC . If the production of reso- 
nances^ and long-range correlations are well described by the simulation, these should 
cancel in the double ratio and only BEC should remain. The agreement between the 
simulation and data was checked and is good for both the unlike-charge and like-charge 
distributions. The latter show significant deviations only in the low Q region, where dis- 
tortions due to BEC are expected in the data. Thus, for the double ratio, a simple fit 
ansatz can be used: 



C'(Q)=N(l + UQ)\e-Q 2R2 ). (16) 

As in section 3.2, the double ratios for the three event selections can be described by 
superpositions of the correlations for the different pion classes. Eqs. [L2| - [n] are also valid 
for the double ratios. It can be shown that the relative probabilities P are given by the 
expressions given in Table |2], except that, in this case, the number of like-charge pairs 
N±± have to be used instead of the number of unlike-charge ones iV + _ as was the case 
in section 3.2. The relative probabilities are determined from a Monte Carlo simulation 
without BEC. 

In a simultaneous fit to the three double ratios C'(Q) the BEC for the three pion 
classes C fsame (Q), C fdlS (Q) and C' Z *(Q) are determined. A common source radius for all 
pion classes is assumed and the parameters A and R in the correlation function C' Z *(Q) 
are adjusted for differences in multiplicity and topology as in section 3.2. The seven 
parameters used in the fit are given in Table [5]. The fit is made in the full range 0.0 < 
Q < 2.0 GeV/c 2 . 

The fit describes the distributions well, with a % 2 /d.o.f. of 72.8/67. The results of the 
fit are given in Table [F[ They are fully compatible with the results of section 3.2. The 
systematic errors are obtained in a similar way as before, with the relevant individual 
contributions given in Table ||. This method has the advantage that the long-range 
correlations do not have to be determined in the fit. On the other hand, this method 
relies more on Monte- Carlo input. 



4 As for the main analysis, the resonance cross-sections in Jetset are adjusted to the measured rates 
at LEP energies. 



Parameter 


same W diff W (Z u /7)* 


R (fm) 
A 
N 


1.11 ±0.13 ±0.21 
0.65 ±0.21 ±0.09 0.50 ±0.78 ±0.14 0.42 ± 0.09 ± 0.05 
0.99 ±0.01 ±0.03 1.00 ±0.01 ±0.00 0.99 ± 0.01 ± 0.02 



Table 7: Result of the simultaneous fit using the double ratio C'(Q) The first error 
corresponds to the statistical uncertainty the second one to systematics. 





R (fm) 


^same 


A ditf 


\ z 


X 2 /d.o.f. 


Reference 


1.10 ±0.11 


0.64 ±0.20 


0.50 ±0.72 


0.42 ±0.09 


72.8/68 


Variation 


6R (fm) 


<5A same 


5A diti 


5X' Ai 




1 


±0.11 


±0.07 


±0.03 


±0.09 




2 


-0.14 


-0.06 


< 0.09 


0.03 


89.6/67 


3 


-0.12 


< 0.01 


±0.06 


±0.03 


91.5/67 


7 


< 0.01 


< 0.01 


±0.01 


< 0.01 


72.8/67 


Total 


0.21 


0.09 


0.14 


0.05 





Table 8: The effect of the systematic variations studied (discussed in Sect. |3.3| ) on the 
variables R, A same , A dlff and A z * from the double ratio. The last column shows the quality 
of the corresponding fit. 

4 Discussion and Summary 



We have analysed the data obtained by the OPAL detector at e + e~ center-of-mass energies 
of 172 and 183 GeV to study BEC between pions in three different physical processes: fully 
hadronic events W + W~ — > qqqq, semileptonic events W + W~ — > qqt&i, and non-radiative 
(Z°/j)* events. The analysis assumes equal source size R for these processes. BEC are 
observed each of these processes. The chaoticity parameter A for the semileptonic process 
W + W~ — ► qqtDi is larger than for the processes W + W~ — > qqqq and (Z°/j)* — > qq, but 
still consistent within the errors. The long-range correlation parameters are consistent 
within their errors. Furthermore, BEC between pions from the same W boson and different 
W bosons have been studied. The result for pions from the same W boson is consistent 
with those for pions from non-radiative (Z /^)* — > qq events. At the current level of 
precision it is not established if BEC between pions from different W bosons exists or not. 
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